Effects of site dilution on the magnetic properties of geometrically frustrated antiferromagnets 
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Abstract 

The effect of site dilution by non magnetic impurities on the susceptibility of geometrically frustrated antiferromagnets (kagome and py- 
rochlore lattices) is discussed in the framework of the Generalized Constant Coupling model, for both classical and quantum Heisenberg spins. 
For the classical diluted pyrochlore lattice, excellent agreement is found when compared with Monte Carlo data. Results for the quantum case 
are also presented and discussed. 

PACS numbers: 75.10.Hk, 75.10Jm, 75.30.Cr, 75.40.Cx 



Introduction.- During recent years, geometrically frus- 
trated antiferromagnets (GFAF) have emerged as a new 
class of n^ftrfi^pinrft^ exhibit novel phases at low 
temperaturesiBBcBBQffi. The elementary unit in these sys- 
tems is based on the triangle, which makes it impossible to 
satisfy all the antiferromagnetic interactions at the same time, 
leading to frustration. Examples of GFAF are the pyrochlore 
and kagome lattices. Experimentally, it is found that the sus- 
ceptibility in materials belonging to this class exhibits a high 
temperature phase in which it follows the Curie-Weiss (CW) 
law. Below the Curie-Weiss temperature, geometrical frus- 
tration inhibits the formation of a long range ordered (LRO) 
state, and the system remains paramagnetic, even though there 
are strong correlations between units. This phase is univer- 
sally present in these systems, and it is called the coopera- 
tive paramagnetic phase. Finally, at a certain temperature, TV, 
which depends on the particular material and, usually, is well 
below the CW temperature, there appear non universal phases: 
some of the systems remain paramagnetic down to the lowest 
temperature reached experimentallyo, some of them exhibit 
non collinear ojjdsred statesQ, or even some of them form a 
spin glass stateBQa, even though the amount of disorder in the 
structure is very small. 

In the simplest theoretical description of GFAF, the spins on 
the lattice are regarded as Heisenberg spins with only nearest 
neighbor (NN) interactions. In this picture, it is predicted that 
the non trivial degeneracy of the ground state inhibits the for- 
mation of a LRO state, and-the system remains paramagnetic 
down to zero temperaturelij. However, due to the presence 
of frustration, the NN exchange does not fix an energy scale 
on the problem, and any small perturbation can break the non 
trivial degeneracy of the ground state and lead to some kind 
of ordered state. Therefore, it is especially important to incor- 
porate these possible perturbations in any model that tries to 
explain the low temperature phases of these systems. Exam- 
ples of perturbations present in real systems are next nearest 
neighbor (NNN) interactions, small anisotropies, long range 
dipole-dipole interaction, or dilution by non magnetic impu- 
rities, to cite some. 

Of course, before trying to understand the non universal be- 
havior observed below XV, it is very important to fully explore 
models that quantitatively describe the universal cooperative 



Figure 1 : Comparison between the inverse susceptibility calculated 
in the GCC model for the diluted pyrochlore lattice for different con- 
centrations of non magnetic impurities x and MC data from Ref. |ll[ 



small clusters can provide a very accurate description of the 
universal cooperative paramagnetic regime in these systems, 
when compared with both Monte Carlo (MC) and experimen- 
tal data. The purpose of the present work is to study what 
are the effects of site dilution by non magnetic impurities of 
the otherwise perfect pyrochlore and kagome lattices on the 
magnetic properties in this same regime. This problem is rel- 
evant for establishing comparisons with experimental data as, 
obviously, any real material contains a certain amount of dilu- 
tion. On more phenomenological grounds, this problem was 
discussed in Refs. [l5| and 11 

The model.- The central idea of the models in 
Refs. p^,pj|Jl4| is based on the experimental and numer- 



paramagnetic phase. In a series of recent papersELjO, the 
present authors have shown how very simple models based on 



ical observation that the correlations in these systems are 
always short ranged. Therefore, we can start by studying the 
properties of a small cluster and, later, add the interactions 
with neighboring spins outside the cluster in an approximate 
way. In the simplest approach, these interactions can be mod- 
eled as a mean field. A more refined approximation consists 
on treating these interactions as created by an effective field 
which is fixed by a self consistency condition, tbsi-sp called 
generalized constant coupling (GCC) approachOliil. This 
self consistency condition is constructed in the following 
way: we consider the magnetization per spin of a unit with 
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p spins (a triangle for the kagome lattice or a tetrahedron 
for the pyrochlore lattice, for example), in the presence of 
the internal field created by the neighboring p — 1 spins 
outside the unit, and the magnetization of an isolated spin 
in the presence of the internal field created by the 2(p — 1) 
neighboring spins, and we equate both quantities, obtaining 
in this way an equation for the internal field 



magnetization for different types of units 

(1 - x) m cw ([H + 2(p - 1)(1 - x)H'] /T) = 

- J2 P2(x) m q {[H + (p - 1)(1 - x)H'] /T) , (6) 

9=1 
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In this expression, Hq is the applied magnetic field and H 1 is 
the internal field; m p represents the magnetization of a unit 
with p spins and m c the corresponding magnetization for a 
single spin in the presence of the internal field created by its 
2(p — 1) NN. Expressions for these quantities can be found 
in references [jj] and |l4| for the classical and quantum cases, 
respectively. Of course, this equation can only be solved nu- 
merically in the general case but, in the paramagnetic limit, 
[Hq + 2(p — 1)H'] /T <C 1, we can expand it up to first order 
in the internal field in terms of the susceptibility of an isolated 
unit and obtain an analytical expression for the susceptibility 
per spin in the presence of the internal field 
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where 



s p (T) = —-tlnZ p (T)-l, 



(2) 



(3) 



with T a dimensionless temperature defined by T j J, where 
J = Jsq in the classical limit (J is the positive antiferromag- 
netic exchange coupling and s the length of the spin) and 
J = Jso( s o + 1) in the quantum case. Z p is the partition 
function of an isolated unit with p spins. Expressions for this 
quantity in the classical limit can be found in Ref. Ill], In the 
quantum limit, this quantity is given byO 
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g{S)(2S + l)e ^oCo+i)? . 



(4) 



where S represents the total spin of the unit and g(S) the cor- 
responding degeneracy. 

Let us now consider the effect of substituting some of the 
magnetic ions in the lattice by non-magnetic ions. We will as- 
sume that the distribution of non-magnetic impurities is com- 
pletely random, so the number of units with q = 1,2,... mag- 
netic ions for a lattice formed by units with p ions, far a con- 
centration of non-magnetic impurities x, is given byliil 



(5) 



The self-consistency condition that determines the internal 
field H' is given, in this case, in terms of the averaged 



where the (1 — x) factor in front of the CW magnetization 
takes into account the reduction of the total number of mag- 
netic ions upon dilution from N to N (1 — x), whereas the 
prefactor in the argument of both magnetizations stands for 
the reduction in the number of NN. 
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Figure 2: Classical susceptibilities for the diluted (a) kagome and 
(b) pyrochlore lattices. The inset shows the same quantity at very 
low temperatures. The curve labeled as "Curie-Weiss" is given by 
1/(1 + T) in this temperature scale. 



Equation (g) can be solved for H', in the paramagnetic re- 
gion, along the same lines mentioned above, by introducing 
the function 



e p (T,x) 



2T 2 3 — ~ 

lnZ p (T,x)-l, (7) 
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which reduces to (||) for x = 0, where Z p (T, x) is an averaged 
partition function defined by 



Z p (T,x) = l[Z q (Tf P «(* 



(8) 



9=1 



With these definitions, eq. (||) can be solved in the param- 
agnetic regime to obtain the averaged susceptibility 



X P (T, x) 



l-xl + e p (T) 



3 JT 1 -e p (T) 
The corresponding internal energy per spin is given by 



(9) 



u p (T,x) = J £ p (T,x), 



and the specific heat 



Cp(T,x) = -~e p (T,x) 
ol 



(10) 



(11) 



Discussion.- The pyrochlore and kagome lattices have been 
extensiiielVj-studied by numerical methods in the classical 
limiOLlLKj. Therefore, there is a wealth of available data 
obtained from MC simulations that allow us to estimate the 
accuracy of our model in the classical limit. Actually, this 
comparison was recently carried out by the present authors 
for the non diluted lattices considered in this work, and it was 
found that both the susceptibility and specific heat are essen- 
tially exact for the pyrochlore lattice down to K. In the case 
of the kagome lattice, however, the susceptibility was found to 
be again essentially exact, but the MC specific heat deviates 
slightly as we go to K, due to the order by disorder effect 
which, obviously, is not included in our model. We will not 
repeat such a comparison here, as it can be found in Ref. [l3| 

There are also some results from MC calculations for the 
case with site .dilution but, unfortunately, only for the py- 
rochlore latticeU. They have been compared with the predic- 
tions of our model in the classical limit in Fig. [l| As can be 
seen, the susceptibility predicted by our model is essentially 
exact down to K. It is important to notice that the curves in 
that figure do not include any adjustable parameter or global 
scale factor. Of course, it would be extremely interesting to 
test the predictions of the model for the kagome case too. 

In Figs. |2| and |3| we have plotted the classical and quantum 
susceptibilities predicted by our model for different concen- 
trations of non magnetic impurities. The quantum suscepti- 
bilities depicted there are for the kagome spin 5/2 and spin 
1, and pyrochlore spin 3/2. The first one would correspond 
to the iron jarosites compounds and the last one to ZnCr204. 
The second one would correspond to Mo 4+ ions on a kagome 
lattice, and has been included to stress the differences be- 
tween integer and half integer values of sq in our model. In 
these figures, we have rescaled the temperature by the CW 

M Ocw(x) = 2(P 3 1)J (1 - x), and the suscepti- 



temperatureli 



bility by the CW constant Ccw {%) 



_ J(l-x) 
3J 



The advantage 

of this rescaling is that all the curves for different concen- 
trations fall on the same curve at high temperatures, which 




Figure 3: Quantum susceptibilities for the diluted (a) kagome (so = 
|), (b) kagome (so = 1), and (c) pyrochlore (so = §) lattices. The 
inset shows the same quantity at very low temperatures. 



makes it easier to identify the nontrivial effects of dilution at 
low temperatures (that is, those which do not come from a re- 
duction of (1 — x) factor in both the CW temperature and CW 
constant). 

From those figures, we can see that the main effect of di- 
lution in the quantum kagome and pyrochlore systems, and 
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also in the classical kagome lattice is that the intensity of the 
maximum in the susceptibility decreases as the concentration 
of impurities increases and, eventually, this maximum disap- 
pears for a certain amount of non magnetic impurities. In the 
classical pyrochlore lattice, there is no maximum in the sus- 
ceptibility, but only a shoulder, due to the formation of short 
range correlations in the tetrahedral unit. Again, this shoul- 
der disappears for a certain amount of non magnetic impuri- 
ties. Another effect of dilution is related to the upturn of the 
susceptibility at low temperatures. Both the classical non di- 
luted kagome and pyrochlore lattices susceptibilities go to a 
finite value at K. In the quantum case, the pyrochlore sus- 
ceptibility in the pure system falls to zero at K. The kagome 
susceptibility falls to zero for integer values of the individual 
spins or diverges for half integer values of sq, even in the ab- 
sence of dilution. This divergence is related to considering a 
cluster with an odd number of half integer spins, which has a 
ground state with a | total spin momentum. This problem of 
the model and its implications has been discussed in Ref. |l4| 
When we add a certain amount of dilution, there is a diver- 
gence of the susceptibility as we approach T = for all the 
cases considered in this work. This upturn in both the classi- 
cal kagome and pyrochlore lattices, and also in the quantum 
pyrochlore and quantum kagome lattices with integer values 
of so, is due to the appearance of loose spins in the system. 
However, in the quantum kagome with half integer values of 
so, this upturn is present even in the pure system, as com- 
mented above. In the quantum pyrochlore lattice with half 
integer sq, both the loose spins and triangular units contribute 
to this upturn of the susceptibility. 



Even though we will not present a detailed comparison with 
experimental measurements here, it is interesting to note that 
some of the features exhibited by the calculated quantum sus- 
ceptibilities are qualitatively.-similar to those experimentally 
found in pyrochlore-systemsoBQ, iron jarositesffl and also in 
SrCrgpGai2_9pOi9Ej (where these last two systems are con- 
sidered as experimental realizations of the kagome lattice), 
namely, a maximum in the susceptibility for both the py- 
rochlore and kagome lattices. The disappearance of this max- 
imum with increasing dilution and the upturn of the suscep- 
tibility below the maximum for pure samples has been also 
observed in the iron jarosites compounds. However, a lot of 
caution has to be exercised when comparing the present model 
with experimental data for two main reasons: first, the region 
where the maximum of the susceptibility is observed, it is pre- 
cisely the region where spin freezing starts to take place and, 
obviously, this effect is not accounted for in our model. Sec- 
ond, it is not obvious that the distribution of non magnetic 
impurities in real systems is purely random and this fact could 
lead to important quantitative differences when trying to use 
this model to interpret experimental data. 

In any case, we think that the present model could provide 
further insight in the properties of the geometrically frustrated 
antiferromagnets in the cooperative paramagnetic region. 
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